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Chapter  IX 
Burner  Flarr.es 

1.  Features  of  an  Anchored  Flame. 

A  common  enough  observation,  both  inside  and  outside  the  laboratory,  is 
that  a  premix^d  flame  can  be  positioned  stably  at  the  mouth  of  a  tube  through 
which  the  mixture  passes.  Such  a  flame  is  usually  conical  (though  r.ct 
necessarily  so),  in  which  case  it  is  conveniently  divided  into  three  parts: 
the  tip,  the  base  (near  the  rim  of  the  tube)  and  the  bulk  cf  the  flame  in 
between.  The  latter  gives  it  roughly  the  shape  of  a  ccr.e. 

Elementary  considerations  of  the  flame  speed  and  the  nature  of  the 
flow  adequately  explain  the  conical  shape  (see  Fig.  VIII. 2  and  the  accompanying 
discussion).  Simple  hydrodynamic  considerations  provide  salient  features  cf 
the  associated  flow  field,  as  we  shall  see  in  Sec.  2. 

The  nature  of  the  combustion  field  in  the  vicinity  of  the  rim  is  crucial 
in  questions  of  existence  and  stability  of  the  flame.  Gas  speeds  near  the  tube 
wall  are  small,  because  of  viscous  effects,  so  that  if  the  flame  could  penetrate 
there  it  would  be  able  to  propagate  against  the  flow,  travelling  down  the  tube 
in  a  phenomenon  known  as  flashback.  In  point  of  fact,  the  flame  is  quenched  at 
some  distance  from  the  wall,  so  as  to  prevent  its  reaching  the  low-speed  region, 
the  fundamental  quenching  mechanism  (for  a  stationary  flame)  being  heat  loss  by 
conduction  to  the  tube.  Such  quenching  enables  unburnt  gas  to  escape  betveeen 
the  flame  and  the  wall  through  what  is  known  as  the  dead  space.  A  mathematical 
description  of  this  phenomenon  is  presented  in  Sec.  5. 

The  heat  loss  to  the  tube  also  plays  an  important  role  in  the  global 
stability  of  the  combustion  field.  If  the  flame  is  disturbed  so  that  its 
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speed  increases  it  will  move  closer  to  the  rim,  threatening  flashback.  But 
heat  loss  will  then  increase,  causing  a  decrease  in  the  temperature  of  the 
flame  and  hence  its  speed,  so  that  it  will  be  swept  away  again.  Similarly, 
if  there  is  a  decrease  in  flame  speed  so  that  the  flame  is  swept  away  from 
the  rim,  threatening  so-called  blowoff,  heat  conduction  to  the  wall  decreases, 
the  speed  increases  and  the  flame  moves  back  again.  Thus  heat  loss  acts  as  a 
restoring  force. 

The  third  part  of  a  burner  flame,  the  tip,  is  similar  in  one  respect  tc 
the  base  in  that  here  also  the  flame  approaches  a  boundary,  albeit  an  adiabatic 
one,  namely  the  centerline.  Experiment  shows  that  the  tip  may  be  closed,  the 
flame  cutting  the  centerline  at  right  angles,  or  open,  the  flame  being 
quenched  (though  not  by  heat  loss)  at  a  finite  distance  from  the  centerline  so 
as  to  leave  a  dead  space  or  hole  through  which  the  mixture  escapes  unburnt 
(Fig.  l).  Which  of  these  two  possibilities  occurs  depends  on  the  composition 
of  the  mixture  and  flow  conditions.  Flame  tips  are  treated  in  Secs.  3  &  k. 


The  hydrodynamic  Jump  conditions  (VIII. 1,2, 3)  are  satisfied  by  the  undisturbed 
flow  (l)  ahead  of  the  flame  and 


(2)  v  *  (U,  ±eU(Pl/p2  -  1)) 


behind,  where  the  sign  is  opposite  to  that  of  the  slope  of  the  flame.  There  is 
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simply  a  refraction  of  the  streamlines,  and  we  may  expect  that  to  be  the 
primary  characteristic  of  the  flow  field  in  the  neighborhood  of  any  point  on 
the  conical  part  of  the  burner  flame.  A  more  detailed  description  requires 
a  hydrodynamic  analysis,  which  we  shall  pursue  some  distance  assuming  that  the 
flame  speed  is  constant.  Such  an  assumption  is  correct  if  the  burner  is 
sufficiently  large;  in  the  context  of  slowly  varying  flame  the  voluminal  stretch 
in  equation  ( VIII. 43)  -is  then  small  so  that  the  burning  rate  =  1. 

We  shall  first  look  for  a  plane  flame  tip  in  which  the  uniform  flows  (l) 
and  (2)  hold  far  upstream;  elsewhere  it  is  necessary  to  solve  Euler's  equations 
on  each  side  of  the  constant-speed  flame.  A  Poisseuille  distribution  would  be 
more  realistic  for  the  upstream  flow  ahead  of  the  flame,  but  the  analysis  is 
much  easier  for  a  uniform  flow  while  still  uncovering  the  essential  characterist 
Although  we  are  not  necessarily  dealing  with  a  slowly  varying  flame,  the 
notation  5,  n  of  Sec.  VIII. 3  will  be  used  to  emphasize  that  the  distance  unit 
is  large  compared  to  the  flame  thickness.  The  £ -axis  is  now  taken  along  the 
symmetry  line  (Fig.  2);  velocity  components  are  u,  v.  Since  the  flame  slopes 
gently  (in  the  main),  we  may  write 

(3)  n  =  F(a)  with  a  =  e  C 

for  its  upper  surface.  In  front  u,  p  and  3/3n  are  0(l)  while  v  and 
3/3£  are  0(  e).  The  momentum  equation  therefore  shows  that  p  is  a  function 
of  o  alone ,  to  order  e  ,  from  which 

(4)  u  *  f(o)  +  0(e  2) 


follows  on  applying  Bernoulli's  equation;  continuously  then  requires 
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(5)  v  =  -  enf ' (a)  +  0 ( e2). 

The  condition  for  the  flame  to  be  at  rest  is 

(6)  fF  =  -Uo  +  0(e2) 

0 

when  a  is  measured  from  the  tip. 

The  appropriate  coordinates  in  the  burnt  gas  are  o  and 

(T)  v  =  en- 

At  the  flame  we  have 

(8)  u2  =  f  +  0(e2),  v2  =  e(Pl U/p2  +  fF' )  +  C(£2) 

according  to  the  jump  conditions  while,  on  the  centerline,  symmetry  requires 

(9)  v  =  0  for  a  >  0. 

Since  the  flow  is  a  perturbation  of  the  uniform  one  upstream,  we  may  write 

(10)  u  =  U  +  eu'  +  0(e2),  v  =  ev'  +  C (£2), 
which  immediately  implies 

(11)  f  =  U  +  0(e)  and  hence  F  =  -o  +  C(£). 

The  flame  is  wedge  shaped  and  the” boundary  conditions  (3)  yield 

(12)  n'2  =  0,  v'2  =  U(p1/p2  -  l) 

if  the  flow  upstream  is  to  be  unperturbed. 


Since  the  vorticity  is  constant  on  the  undisturbed  streumlines,  the 

flow  downstream  is  the  superposition  of  a  potential  motion,  satisfying  the 

boundary  conditions  (9)  and  (12b),  and  a  shear  motion  (with  v'  =  C)  induced 

by  the  boundary  conditions  (12a).  Thus 

U(p  /p  -  1)  U(p  /p  -  1) 

(13)  u'  =  — - - -  i.n(v/er),  v'  =  - $ 

If  7T 

where  r,4>  are  polar  coordinates  corresponding  to  a,v  •  The  flow  field  is 
sketched  in  Fig.  2.  Note  the  singularity  in  u'  for  <j>  =  0,  which  presumably 
implies  that  the  flow  differs  significantly  from  uniform  there.  The  nature 
of  this  wake  has  not  been  investigated,  so  that  the  description  is  incomplete 
at  the  present  time. 

A  similar  treatment  is  possible  for  axisymmetric  flames;  only  the  changes 
in  the  previous  formulas  need  be  mentioned.  With  n  denoting  distance  from 
the  axis  of  symmetry,  the  corresponding  velocity  component  is 

(ll+)  v  =  -  ~  epf'(a)  +  0(e  ), 

so  that  the  flame  is  at  rest  if 

(15)  fF'  +  -f'F  =  -U  +  0(e2). 

Nevertheless  the  results  (ll)  still  hold,  as  do  the  boundary  conditions  (9) 
and  (12).  We  now  find 

(16)  u'  =  0,  v'  =  j  e  U(p1/p2  -  l)tan  4>/2, 

showing  that  the  disturbance  is  o(e)  everywhere  except  near  the  cone 
^  s  ir-e  on  which  the  flame  lies.  That  is  reflected  in  Fig.  3,  where  the 
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streamlines  straighten  out  much  faster  than  in  Fig.  2;  there  is  no  wake. 

Note  that  the  deflection  of  the  streamlines  at  the  flame  is  the  same  in  the 
two  cases,  namely  elKp^/p^  -  l). 

The  streamline  tracing  of  actual  flames  (Fig.  l)  differ  in  two  respects. 
The  flow  ahead  of  the  flame  is  not  uniform  but  is  essentially  Poisseuille; 
and  finite  values  of  e  give  rise  to  a  pressure  gradient  which  is  associated 
with  the  divergence  of  the  streamlines  there. 


3.  Slowly  Varying  Flame  Tips. 

The  question  that  we  examine  here  is  whether  the  equation  (’/III. 32) 
governing  slowly  varying  flames  admits  solutions  corresponding  to  the  tips 
(Sivashinsky  197^,  1975)-  The  difficulty  cf  the  hydrodynamic  problem,  ever, 
when  the  burning  rate  is  prescribed  (Sec.  2)  ,  leads  us  to  sever  the  coupling 
with  Euler's  equations  by  setting 


(17)  y1  =  u(i,o,o), 


in  accordance  with  the  constant-density  approximation  (Sec.  1.5). 
If  iji  is  the  inclination  of  the  flame,  then 

(l8)  =  p^U  sin  ip,  yj^  =  U  cosip 

so  that,  in  the  plane  case  (Fig.  4), 


(19) 


Xu  'hMn  =  pl^cosMf’  Y±.'Xli 


and  the  governing  equation  becomes 


(20) 


^  =  k  sin\fcn(-Sv"  , 
ds  r  sin  a 
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where 

(21)  k  =  -2T?pfu/b  and  a  =  cocec  ^p. U. 

*  1  .1 

The  implicit  requirement  p^U  >  1  simply  reflects  the  fact  that  a  stationary 
flame  can  only  exist  in  a  flow  that  moves  faster  than  the  adiabatic  flame  speed. 
For  comparison  with  the  axially  symmetric  case  to  follow,  it  is  more  convenient 
to  use  n  as  independent  variable,  in  terms  of  which  the  equation  becomes 

(22)  k  sin2«Hn(4^). 

dp  sm  ct 


Apart  from  =  a  and  rr  -  a  ,  only  three  integral  curves  are  shown  in 
Fig.  4a.  The  remainder  can  be  found  by  horizontal  translation.  These  which 
lead  to  shapes  with  $  -*■  a  or  it  -  a  as  £  ,  as  required,  are  shewn  in 

Figs.  4b, c.  In  Fig.  4b  the  curves  clji  and  (5)  represent  a  family  filling  in  the 
wedge  between  the  straight  lines  (2}  and  (4)  .  In  Fig.  4c,  only  these  lines 
are  obtained. 

For  g£>  1  (i.e.  b  <  0  and  k  >  0)  a  smooth  tip  is  possible,  namely  (]Q  , 
and  presumably  that  should  be  chosen  over  the  whole  family  of  cusped  tips  which 
are  available.  For  <  1  (i.e.  b  >  0  and  k  <  0) ,  no  smooth  tip  is  possible; 
only  the  wedge-shaped  tip  formed  by  (2".  and  is  offered  by  the  theory. 

Turning  to  the  axisymmetric  case,  we  find  that  equations  (l8)  and  (19a) 
still  hold  hut  that  (19b)  is  replaced  by 


(23) 


V'll  =  n  h  (nCOS  *} 


because  of  the  radial  (p)  divergence.  Hence  the  governing  equation  now  becomes 


(24) 


=  k  sin2^n(4ni), 

dp  n  sina 
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vith  k  and  ot  defined  as  before,  where  n  ranges  over  positive  values 
only. 

Because  of  the  new  term  the  integral  curves,  shown  in  Fig.  5a,  no  longer 
have  the  translation  property,  though  they  have  the  same  shape  as  in  the  plane 
case  as  n-*”  •  Cnly  those  asymptoting  ij  =  i  -  a  as  are  acceptable,  which 

leads  to  the  open  (including  re-entrant)  possibilities  in  Fig.  1+b,  as  well  as 
one  closed,  but  only  to  the  single  open  possibility  in  Fig.  Uc.  Representative 
curves  are  shown,  as  marked;  all  open  curves  approach  the  centerline  as 

5  ■*  ±  <». 

For  J_>  1  (i.e.  b  <  0  and  k  >  0)  a  closed  tip  is  possible,  namely  Q  , 
and  presumably  that  occurs  rather  than  either  a  re-entrant  or  open  one.  ( Cf . 
later  remarks  concerning  flame  temperature.).  For  <1  (i.e.  b  >  C  and  k  <  l) 
only  one  tip  is  available  and  that  is  open.  In  short,  the  theory  of  slowly 
varying  flames  predicts  closed  tips  for  .'  >  1  ar.d  open  tips  for  ^_<  1.  It  is 
noteworthy  that  open  tips  are  found  experimentally  (lewis  &  von  Elbe  19^1,  p.  21 
for  both  lean  hydrogen  and  rich  heavy  hydrocarbon  mixtures.  In  each  case 
the  deficient  component  diffuses  much  more  readily  than  the  component  in  excess 
and  indeed  calculations  based  on  values  given  by  Kanury  (1977,  pp.  385  et  sec.) 
show  that  <  =  0.23  and  0.8l,  respectively,  for  mixtures  with  air. 

The  flame  temperature  is  readily  deduced  from  the  formula  ( VIII. 30 )  as 

(25)  H.,  +  20  in  M  =  H  +20  ^Hfin(sinii)/sina) . 

1  1  n  1  1 

In  all  cases  the  adiabatic  flame  temperature  is  approached,  of  course, 

as  ip-*-  Tt-a  at  large  distances.  When  7)  >  1 ,  the  temperature  tends  to  a 


maximum 
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+  20  cosec  a 

at  the  very  tip,  irrespective  of  whether  the  flame  is  plane  or  axisymmetric . 

On  the  other  hand,  when  ,1  <  1  the  temperature  does  not  change  along  a  plane 
tip  hut  decreases  like 

(2?)  H  +  2e_1f{^n  * 

for  an  axisymmetric  tip,  giving  an  unbounded  perturbation.  (The  same  result 
holds  for  the  open  tips  in  Fig.  5b,  so  that  the  closed  tip  selected  is  the 
only  one  without  this  singularity. )  Such  unboundedness  means  that  the 
temperature  differs  by  more  than  0(8  from  the  adiabatic  value;  so  the 
asymptotic  structure  is  not  uniformly  valid  and  a  different  kind  of  analysis 
is  required.  While  such  an  analysis  has  not  yet  been  attempted,  the  drop 
in  temperature  suggests  that  the  flame  is  extinguished;  so  that  the  profile  cf 
Fig.  kc  should  be  truncated  at  some  large  value  of  £  ,  giving  a  flame  core  in 
accord  with  experimentally  observed  open  tips.  For  the  plane  tip,  non¬ 
uniformity  manifests  itself  geometrically  rather  than  through  a  singularity 
in  the  flame  temperature:  the  pointed  tip  violates  the  assumption  of  slow 
variation  on  which  the  analysis  is  based. 

The  dramatic  change  in  the  nature  of  the  solution  as  passes  through  1, 
together  with  the  breakdown  of  the  analysis  of  1  =  O(0-^),  suggests  that 
flame  tips  should  be  investigated  when  the  diffusion  of  heat  is  almost  the  same 
as  that  of  the  reactant,  i.e.  under  the  near-equidiffusional  assumption.  That 
has  been  done  by  Buckmaster  (1979b),  whose  work  is  the  subject  of  the  next 


section. 


4 .  Near-Fquidif fusional  Flame  Tips. 


The  general  formulation  of  near-equidiffusional  flames  as  elliptic  free¬ 
boundary  problems  is  given  in  Sec.  VIII. 5.  Moreover,  if  the  constant-density 
approximation  is  invoked  to  uncouple  the  fluid  mechanics  and  the  flame  is 
immersed  in  a  uniform  flow  whose  speed  is  much  greater  than  the  flame  speed, 
then  these  problems  reduce  to  ones  of  Stefan  type  (as  is  shown  in  Sec.  VIII.  6). 

To  describe  flame  tips,  take  the  x~axis  along  the  centerline.  For  upstrea 
the  flame  is  either  wedge  shaped  or  conical  with  half-angle  it  /U,  corresponding 
to  the  adiabatic  value  1  for  the  flame  speed.  Only  y  >  0  is  considered, 
because  of  symmetry  in  the  plane  case  and  by  definition  in  the  axially  symmetri 
Take 

(28)  y  *  F(x) 

as  the  locus  of  the  flame,  choosing  the  origin  of  x  sufficiently  far  upstream 
for  F(0)  to  be  large.  Hear  x  =  0  the  combustion  field  is  that  of  a  plane 
flame,  so  that  equations  (ill. 13,1*+)  hold  with  ^  =  1  +  X/0  and  x  measured 
vertically  down,  i.e.  on  x  =  0  we  have 

(29)  T  =  T1  +  Y1ey-F^0^  h  =  XYJFvO)  -  y]ey_F^  for  0  <  y  <  F(0), 

(30)  T  =  T*  =  Hx,  h  =  0  for  y  >  F(o), 

which  will  be  taken  as  initial  conditions.  Stricly  speaking,  these  conditions 
are  valid  only  in  the  limit  F(0)  -*■  ®  but,  because  the  exponential  can  be  made 
arbitrarily  small  near  y  =  0,  any  desired  accuracy  can  be  obtained  by  taking 
F(0)  sufficiently  large. 
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The  equations  to  be  solved  are 


(31)  3T / 3 X=  (L  (t)  ,  3h/3x  =  *Lv(h  +  XT) 


(32)  T  =  T#  =  H  ,  3h/3x  =  (L^h)  for 


(cf.  Sec.  III. 6),  where 


(33) 

v  yv  3y  sy 


for  0  <  y  <  F(x), 
y  >  F(X) 


and  v  =  0  or  1  according  as  the  flame  is  plane  or  ax i symmetric.  The 
accompanying  jump  conditions  are 


(3M  [3h/3y]  =  -X[3T/ay]  =  XY1exp(hJt/2T^)  at  y  =  F(x), 

T  and  h  being  continuous  there.  To  these  we  must  add  the  symmetry  conditio 


(35)  3T/3y  =  3h/3y  =0  at  y  =  0 


and,  finally, 

(36)  h  -*■  0  as  y  , 

consistent  with  the  formula  (30). 

At  first  glance  it  seems  questionable  that  the  present  formulation,  based 
as  it  is  on  the  assumption  that  the  gas  speed  is  much  greater  than  the  flame 
speed,  could  provide  an  adequate  description  of  flame  tips.  Indeed  at  a 
smooth  closed  tip  the  two  speeds  are  equal  on  the  centerline.  But  there  is 
no  such  difficulty  for  open  tips;  if  the  flame  speed  is  small  in  the  far 
field,  we  may  expect  it  to  remain  small  everywhere  so  that  uniformly  valid 
results  are  obtained.  Our  analysis  could  therefore  lead  to  open-tip  solutions 


only,  whatever  the  parameter  values.  In  fact  it  also  generates  closed- 
tip  solution's,  albeit  ones  that  are  not  valid  near  the  centerline.  Since  the 
unit  of  length  is  the  preheat  thickness  (at  x=  0) ,  the  region  of  non-uniformity 
is  extremely  small. 

The  analysis  is  completed  by  numerical  integration  which  for  such  parabolic 
systems  involves  marching  downstream,  i.e.  the  direction  of  the  time-like 
variable  x  •  Both  plane  and  axisymmetric  flames  have  been  computed  by 
Buckmaster  (1979),  for  T^  =  0.2,  =  1.0  and  several  values  of  X.  He  also 

computed  the  limit  solution  as  X-v  -°°  in  each  case,  to  which  we  shall  come 
shortly. 

Results  for  plane  flames  are  shown  in  Fig.  6.  Note  that  there  is  no 
appreciable  curvature  until  F  has  decreased  to  about  3,  at  least  for  X  <10. 

Then  the  cold  fringe  of  the  pre-heat  zone  reaches  the  centerline  and  inter¬ 
action  with  the  conditions  (35)  begins,  the  effect  of  the  interaction  depending 
on  X  .  For  non-negative  values  the  speed  of  the  flame  increases  monotonically , 

bending  it  around  towards  the  centerline,  which  it  intersects  at  right  angles. 

1  /2 

Apparently  the  tip  is  smooth  with  F  proportional  to  (xQ  -x)  7  near  the 
intersection  point  xQ>  a  conclusion  that  is  supported  by  a  similarity  solution 
there  for  X  =  0  (Buckmaster  1979a).  These  values  of  X  give  1,  when  smooth 

closed  tips  were  also  found  for  slowly  varying  flames.  Moreover,  as  X  increases 
above  10  the  flame  begins  to  curve  at  progressively  larger  values  of  F,  as  if 
it  were  trying  to  get  onto  the  scale  of  the  slowly  varying  flame. 

For  small  enough  negative  values  of  X  the  flame  still  closes  smoothly 
although  its  speed  initially  decreases,  resulting  in  a  somewhat  extruded  tip. 

For  larger  negative  values  the  flame  speed  decreases  to  zero  at  the  point 


labelled  Q  where  the  curve  has  a  horizontal  tangent.  The  curve  for 
X  =  -50  typifies  this  behavior,  which  is  shown  in  extreme  form  for  A= 


—  03. 


The  limiting  curve  as  A-*  -<*>  can  be  obtained  as  follows.  Assume  that  h 
is  0(  A)  for  0  <  y  <  F( y)  but  remain  0(l)  for  y  >  F(  y) ,  as  is  suggested  by 
the  initial  conditions  (29).  Then  the  function 

(37)  h  =  h/X 

is  identically  zero  for  y  >  F(y),  i.e.  outside  the  tip,  and  we  need  only 
solve  the  simpler  problem 

(38)  =  lv(t)>  =  L/h+T)  for  0  <  y  c  F(x)’ 

(39)  T  =  +  Y1ey_F(0)>  h  =  Y1[F(0)-y]ey'?(°)  at  x  =  0  for  0  <  y  <  «(0), 

(1*0)  T  =  H1  ,  h  =  0,  3(T  +  h)/ay  =0  at  y  =  F(y), 

(1*1)  3T/3y  =  3h/sy  =  0  at  y  =  0 

inside  the  tip. 

The  perturbation  temperature  h#  tends  to  follow  the  changes  in  flame 
speed.  It  has  a  finite  negative  value  at  Q,  showing  once  more  that  curved 
flames  cannot  necessarily  be  treated  as  locally  plane  (when  h#  would  tend 
to  -<*>).  The  numerical  integration  can  be  continued  downstream  from  5. 
without  difficulty,  the  profile  first  curving  away  from  the  centerline  and 
then  back  to  it  to  form  a  closed  tip.  There  is  no  experimental  evidence 
for  such  a  bulbous  tip  and  indeed  having  burnt  mixture  ahead  of  the  flame,  as 
would  be  the  case  Just  past  Q,  would  certainly  need  Justification.  Our 
first-order  analysis  of  the  flame  sheet  does  not  distinguish  between  the 


Jk 
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two  sides:  there  are  no  convection  terns  in  equations  ('/III. 55).  It  is 
conceivable  that  at  the  next  order,  which  does  take  account  of  convection, 
there  is  no  acceptable  structure.  Such  would  be  a  sound  nathematical  reason 
for  rejecting  the  solution  beyond  \  =  x0* 

Buckmaster  proposes  to  ternir.ate  the  solutions  at  so  that  Q,  is 

thereby  identified  with  quenching.  Icwnstream  there  is  no  chemical  reaction, 
the  unburnt  gas  below  Q  simply  mixing  with  the  burnt  gas  above.  That  is, 
the  values  of  T  and  h  on  x  =  Xr  are  use<i  as  initial  conditions  for  a 
smooth  solution  of  the  parabolic  equations  (31),  which  now  hold  for  all 
values  of  y.  Such  a  reacticnless  solution  is  consistent  in  that  T 
immediately  drops  below  H^,  the  temperature  at  which  reaction  occurs:  while 
the  initial  temperature  for  y  >  F(xq)  is  elual  to  E^,  that  for  0  <  y  <  ~(xq 
is  less.  Agreement  with  slowly  varying  flames  for  7l<  1  can  be  obtained  by 
supposing  that  the  wedge  shape  found  there  is  actually  open  at  the  vertex. 

The  axisymmetric  results  (Fig.  7)  are  qualitatively  similar,  tut  there 
is  one  sharp  quanititative  difference.  Azimuthal  curvature  affects  the  flame 
speed  (through  the  terms  in  y  ^3/3y)  long  before  the  pre-heat  zone  reaches 
the  centerline.  As  a  consequence,  the  flame  starts  to  curve  much  further  from 
the  centerline  (at  about  F  =  *-0)  than  for  the  plane  flame.  (Figs.  6  &  7 
should  be  compared. )  Likewise  the  dead  space  between  Q  and  the  centerline 
is  more  than  five  times  larger.  In  short,  open  plane  tips  are  much  narrower 
than  their  axisymmetric  counterparts,  a  result  also  .found  for  slowly  varying 
flames  with  1  if  the  wedge  is  supposed  to  be  open  on  the  present  scale. 

A  similar  comparison  holds  for^>  1,  so  that  the  present  analysis  provides  a 
smooth  transition  from  the  general  characteristics  of  slowly  varying  flames 
for  £<  1  to  those  for  £  >  1. 
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A  word  of  caution  is  needed.  The  transition  from  a  closed  to  an  open  tip 
occurs  for  large  negative  values  of  between  -20  and  -50  for  a  plane  flame 
and  -50  and  -ICO  for  an  axisymmetric  one.  Since  A/0  is  assumed  small,  the 
result  will  probably  be  quantitatively  meaningful  only  for  large  values  of  the 
activation  energy  6  ,  values  so  large  as  to  be  seldom  (if  ever)  realized  in 
practice. 

5.  Quenching  by  a  Cold  Surface. 

A  crucial  r.ole  in  the  existence  and  stability  of  a  burner  flame  is  played 
by  quenching  near  its  base,  as  we  have  already  remarked  in  Sec.  1.  However, 
the  situation  near  the  rim  of  a  burner  is  complicated,  involving  not  only 
heat  transfer  but  also  a  multi-dimensional  flow  field  with  mixing  between,  the 
reactants  and  the  cold  ambient  atmosphere.  To  threw  some  light  on  the 
phenomenon  cf  quenching  by  a  cold  wall,  the  present  section  will  examine  a 
simple  model  ( Buckmaster  1979a)  which,  admittedly,  discards  much  that  is 
relevant  to  the  combusion  field  near  a  burner  rim.  Since  heat  transfer  by 
conduction  plays  the  central  role  in  the  phenomenon,  near-equidiffusional 
flames  provide  an  appropriate  framework. 

Consider  a  parallel  flow  of  combustible  gas  moving  over  the  plane  y  =  C 
with  a  velocity  U  .  A  linear  shear  flow  is  the  simplest  representative  of 
realistic  flows  near  the  surface,  but  any  other  could  be  handled  ^ust  as 
easily  (Buckmaster  also  considers  a  uniform  flow.)  Located  in  the  flow  is 
a  stationary  premixed  flame  (Fig.  8),  which  is  influenced  by  the  presence 
of  the  wall  as  a  diffusion  inhibitor  and  heat  sink. 

The  flame  locus  is  again  described  by  equation  (28)  with  the  origin  of 
chosen  to  make  F(o)  large.  Far  upstream,  where  the  flame  is  uninfluenced 
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by  the  wall,  it  still  propagates  with  unit  velocity  but  now  into  a  flow  wj.th 
speed  UF(x).  Consequently  it  follows  the  parabola 

(J*2)  F  =  [F2(0)  -  2x]1/'2, 

though  this  fact  does  not  change  the  initial  conditions  (29,  30)  but  is  solely 
a  check  on  the  subsequent  computations. 

In  focussing  on  -the  effect  of  heat  loss  we  shall,  for  simplicity,  take 
1  =  0.  In  the  limit  U  -*■  ■*>  the  equations  with  which  we  have  to  deal  are 

(1*3)  y3(T,h)/3x  =  L  (T,h)  for  0  <  y  <  F(x), 

(hk)  T  =  T„  =  H  ,  y3h/3x  =  IL  (h)  for  y  >  F(x), 

these  being  the  modifications  of  equations  (31,32)  which  take  account  of  the 
shear  flow.  The  accompanying  jump  conditions  are 

(1*5)  [3T/3y]  =  -Y1exp(hJt/2T»),  [3h/3y]  =0  at  y  =  F(x) 

while  T  and  h  are  continuous  there.  The  condition  (36)  still  holds,  of 
course. 

Finally  we  come  to  conditions  at  the  wall.  The  heat  flux  there  is 
assumed  small  in  order  to  avoid  excessive  variations  in  the  flame  temperature 
(cf.  Sec.  VIII. 2);  specifically 

(U6)  3T/3y  =  k(T  -  T  )/6  at  y  =  0 

to  all  orders,  where  k  is  a  prescribed  positive  constant  and  is,  as 

before,  the  temperature  of  the  fresh  mixture.  Such  a  condition  is  appropriate 
for  a  thin,  poorly  conducting  wall  whose  other  surface  is  maintained  at  the 
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temperature  T  .  Otherwise  it  may  be  considered  an  assumption  which  avoids 
coupling  the  temperature  field  in  the  wall  with  that  in  the  flow  without 
compromising  the  essential  physics.  We  shall  also  take  the  wall  to  be 
chemically  inactive,  i.e. 

(47)  3Y/3y  =0  at  y  =  0 


to  all  orders.  The  boundary  conditions  are  therefore 


(48)  3T/3y  =  0,  3h/3y  =  k(T  -  T.,)  at  y  =  0. 


The  main  simplification  for  1=0  is  that  h  can  be  removed  from 
the  problem.  Thus  y3h/3x  =  IL0(h)  everywhere  and  h  has  no  discontinuities 
at  the  flame  sheet,  so  that  it  can  be  written  explicitly  in  terms  of  its 
normal  derivative  at  the  wall  and  hence,  according  to  the  condition  (48),  the 
surface  temperature.  At  the  flame  sheet  it  becomes 


31/6r(i/3)  /  fyrts.QH  i  f3(x)] 

L  5  *-s 


which  then  determines  the  temperature  gradient  (45)  on  the  cold  side  of  the 


flame  in  terms  of  the  surface  temperature  up  to  that  station.  We  are  left 
with  a  (numerical)  problem  for  T  and  F  alone. 

Results  are  shown  in  Fig.  8  for  various  values  of  the  heat-loss  paramete 

The  flame  sheet  intersects  the  wall  at  right  angles  when  k  is  small,  F 

1/3 

being  proportional  to  (xq  ~x)  near  the  intersection  point  Xq  (as 
Buckmaster  shows  by  a  similarity  solution  there).  Exactly  what  happens  very 


close  to  the  wall  cannot  be  revealed  by  our  analysis,  since  the  assumption  of 
a  large  flow  speed  (which  reduces  the  elliptic  problem  to  a  parabolic  one)  is 
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not  valid  when  y  is  small.  A  reasonable  conclusion  from  the  failure  of 
the  heat  loss  to  generate  a  dead  space  oh- the  scale  of  the  flame  thickness 
is  that  flashback  occurs. 

When  k  is  large  enough  7-2ir/3^^r(l/3)  is  sufficient)  the  flame 

sheet  becomes  horizontal  at  some  point  Q,  which  is  identified  with  quenching 

as  in  Sec.  4.  (Again  the  numerical  solution  can  be  continued  beyond  0  until 

ultimately  the  flame  sheet  intersects  the  wall. )  Failure  to  adopt  such  a 

hypothesis  in  the  present  context  would  lead  to  the  unexpected  conclusion 

that  0(0  heat  loss  is  incapable  of  quenching  a  flame.  We  therefore 

propose  that  the  region  between  Q  and  the  wall  is  a  deadspace  through  which 

unburnt  mixture  can  pass  downstream,  the  generation  of  such  a  gap  being  the 

mechanism  by  which  heat  loss  prevents  flashback.  The  conclusion  is  ir.  sharp 

contrast  with  that  for  a  plane  flame  (Sec.  III.1*),  for  which  extinction  occurs 

-1/2 

when  the  loss  is  sufficient  to  reduce  the  flame  speed  to  e 


j 
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